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Abstract 



We study the long wavelength limit of a spin S Heisenberg antiferromagnetic chain. 
The fermionic Lagrangian obtained corresponds to a perturbed level 25 SU(2) Wess- 
Zumino-Witten model. This effective theory is then mapped into a compact U(l) 
boson interacting with Z2S parafermions. The analysis of this effective theory allows 
us to show that when S is an integer there is a mass gap to all excitations, whereas 
this gap vanishes in the half-odd-integer spin case. This gives a field theory treatment 
of the so-called Haldane's conjecture for arbitrary values of the spin S. 
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Typeset using REVI^X 



In 1983 Haldane conjectured that in the case of integer spin, the spin 5 quantum Heisen- 
berg Antiferromagnetic (HAF) chain has a unique disordered ground state with a finite 
excitation gap, while the same model has no excitation gap when 5 is a half-odd-integer 
|Ij]. Using a mapping to the non-linear o Model, valid in the large 5 limit, the origin of the 
difference has been identified as being due to an extra topological "0 term" in the effective 
cr-model Lagrangian for systems with half- integer 5 ffl. This clearly suggests that the origin 
of this difference is non-perturbative. Although Haldane's predictions were based on large-5 
arguments, it is known that this conjecture is consistent with the Bethe Ansatz exact solu- 
tion which is available for 5 = 1/2 0, and experimental, numerical, and theoretical studies 
almost confirmed its validity for 5 = 1 ||. For higher half-odd- integer spin the ground 
state is either degenerate or has a massless excitation |3| which suggests but does not prove 
critical correlations. In || massless behavior in 5 = 3/2 has been tested numerically (see 
also ||). Recent experimental evidence of the existence of the Haldane gap for 5 = 2 HAF 
chains has been found in the study of the compound MnCl 3 (bipy) [[/]. 

Bosonization techniques have been extensively used in the study of spin chains (see || for 
a review on the subject and references therein). The starting point consists on the mapping 
of the original problem, in terms of spin variables, to the problem of fermionic variables 
with additional constraints (which are necessary for the two systems to be equivalent) and 
then bosonize the resulting system. In ||, the HAF model was formulated as a certain limit 
of the Hubbard model. By using a bosonization technique and a renormalization group 
analysis, they found an effective theory for the low-energy physics of 5 = 1/2 and the limit 
of large 5 spin chains. In [nj another approach was used which relies on the fact that 



the spin 5 HAF chain can be represented as 25 spin 1/2 chains for large ferromagnetic 
coupling. In a recent paper [JIT] the issue of implementing explicitly the above mentioned 



additional constraints within the path-integral approach was taken up, and a reconfirmation 
of Haldane's conjecture was obtained for large values of the spin 5, by using non-abelian 
bosonization techniques. 

In this Letter we study the effective low energy theory in terms of a fermionic coset model 



which corresponds to the level 25 SU (2) Wess-Zumino-Witten (WZW) theory [0,0 follow- 
ing and map it into a coupled system of a compact U(l) boson and Z 2 s-Parafermions 
(PF). The analysis of the phase diagrams of Z<is models together with the knowledge of their 
operator product algebra provide us with the elements to establish the difference between 
integer and half-integer spin chains. 

We would like to stress that our analysis does not rely on a large 5 approximation, but 
is valid for all values of 5. 

Let us briefly review the work of flT|| : The spin 5 HAF model can be written in terms 
of fermionic operators C a i X , with a =|, |; % = 1, 25. a, % and x are spin, color and site 
indices, respectively. The spin 5 operator on one site x is represented by 

c _ /~rt a ° L l 3 n h ^ 

where a a p are the Pauli matrices. In order to correctly represent the spin 5 chain, the 
physical states must satisfy 



J2Ct ix C aix \phys >= 2S\phys 



> 
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J2cL^C ajy \phys>=0, (2) 



ft T a 
ij 



where r a are the SU(2S) generators. The first constraint imposes the condition that allows 
only one spin per site, whereas the second one states that the physical states must be color 
singlets. 

The Heisenberg Hamiltonian 

= E ^.Sy (3) 

<xy> 

can be expressed in terms of the fermionic operators introduced above as 

H = ~ E ClixCajyCpjyCpix + constant. (4) 

^ <xy> 

which has a local SU(2S) x U(l) invariance. This quartic interaction can be rewritten by 
introducing an auxiliary field B as 

H = \ E (Bgcl x c ajy + h.c. + b^b%). (5) 

Z <xy> 

In the mean field approximation B is a constant 2S x 2S matrix, and H can be diagonalized. 
To obtain an effective low energy theory we keep the operators C(k) with k near the Fermi 
surface ±7r/2a. We can then write 

~^C aix = e^ x ^ Rai {x) + e-% x q Lai (x), (6) 



where ^ appears for dimensional reasons and ^ R,La.i{%) are slowly varying on the lattice 
scale. 

We also expand the field B as 



B xy = B e aV *y ~B (1 + aV xy ), (7) 

and define A 1 = \(V xy — V^ y ) and R xy = \{V xy + V^ y ). When substituted back into the 
Hamiltonian, the expansion ([?]) leads to a quadratic integral in R xy which can be performed 
to give 

H = B (-itf m (6 ij d x + Al j )* Rj + 

'%,('\A- + 4)*^) + \ - ^R^Lj) 2 , (8) 

The effective Lagrangian is obtained by introducing a Lagrange multiplier A in the Lie 
algebra of U(2S), together with the use of the identity 

S[^M = lim e - X2 I d2x ( *^ )2 , (9) 



in order to implement the constraints (0) (see 11] for details). The effective Lagrangian 
then reads 
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L = V>fiDJl! - Ai(i^75^f - A 2 (i*i*j) . ( 10 ) 

where _D M = —ia^ + B^, and we have decomposed, for later convenience, the field into 
a 17(1) field a„ and a S77(2£) field 
The Lagrangian can be rewritten as 

L = ^ ia Yi(dt, - ia^ijSap + B^5 aP )^ j(} + 
+ 4(Ax + A 2 ) J R . J l + (Ai + X 2 )3r3l + 

- (Ai - A 2 )(*L a *I, i a*Sy / 9*i0 + fc-c), (11) 

where Jr iL = ^ifa^*^^, jj?,L = **Si,£ia*fl,Lta are S77(2) 2 s and U(l) currents respec- 
tively. 

The first term in the Lagrangian ( |TTD corresponds to the fermionic coset version of the 
level 2S SU(2) WZW theory |14| as was already observed in [11]]. In this context the original 



spin operator corresponds to the fundamental field of the WZW model. The third term can 
be absorbed by a redefinition of the £7(1) gauge field a M . 

Thus, we have to deal with the second and last terms in (TTT|) which can be expressed as 



AC = (A x - \ 2 )(g a 0p a + h.c.) + 4(Ai + A 2 ) f R J L (12) 

where g = ^> Ri $> Li is the spin 1/2 primary field of the SU(2) 2 s WZW theory, & l / 2 \ with 
conformal dimensions h = h = 3/(8(5' + 1)). The first term in ( |T2| ) corresponds to the spin 
1 affine primary $^ with conformal dimensions h = h = 1/(S + 1), so we can write 

AC = -4 (Ai - A 2 ) tr + 4(A X + A 2 )7r.7l (13) 



The 5=1 case is simpler, as has been discussed in 0, [|TTJ, since affine (Kac-Moody) 



selection rules forbid the appearance of the relevant operator $W. We then have an effec- 
tive massless theory in accordance with Haldane's predictions. The second term in ( |i~3"D is 
marginally irrelevant since Ai + A 2 is positive, and gives the well-known logarithmic correc- 
tions to correlators. 

For higher spins, we have to consider the interaction term ([13]) and we also have to 
include all other terms which will be radiatively generated. We then need the operator 
product expansion (OPE) coefficients among the different components of which have 
been computed in jn|. The OPE coefficients are non- vanishing iff the so called "Fusion 
Rules" are non- vanishing. In the level k SU(2) WZW theory they are given by fL6| 



min(j+j',k—j—j') 

$&')_ x $(?') $(n) (u \ 

n=\j-j'\ 

We will use that g0 

ST7(2) fc = Z k ® Z7(l) (15) 

in the sense that the Hilbert spaces of the two theories coincide. We will exploit this 
equivalence to derive an effective low energy action for the spin S HAF chain. Indeed, it 
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was shown in [15| that the primary fields of the SU(2)k WZW theory are related to the 



primaries of the Z^-parafermion theory and the U(l) vertex operators. They are connected 
by the relation 

z) = z) : e ^^m) . (16) 

where the $ fields are the invariant fields of the SU(2)k WZW theory, the fields are the 
Zk PF primaries and ip and (p are the holomorphic and antiholomorphic components of a 
compact massless free boson field. In the same way, the currents are related as 

J+(z) = (2S)^Mz) ■ exp (-^¥>(*)1 : > 



2S 

J° R (z) = (2SY/ 2 dMz) (17) 

where ipi is the first parafermionic field. (A similar relation holds for the left-handed cur- 
rents). 

Using this equivalence we can express the relevant perturbation term fllTf ) as 

,(2) , , (2) . i^ z )-^z)) . 



AC = -4(A X - A 2 ) 0$ + : evm 



+0% : e -vfe(^«) 



+AS{\ 1 + A 2 ) UiV'I = e15S mg, - nt,) : +/i.c.J , (18) 

where we absorbed the derivative part of the U(l) field coming from ( |i~7j ) into a redefinition of 
the constant in front of the unperturbed Lagrangian. The first term corresponds to the first 
"thermal" field of the PF theory, 0q 2 q = e\, with conformal dimensions h = h = 1/(1 + S), 
while the second and third terms correspond to the p = 2 disorder operator in the PF theory, 
02^-2 = ^2 an d its adjoint 0*?2.2 = lA with dimensions d 2 = di = (S — 1)/(2S(S + 1)). 

Including all the operators which are radiatively generated we get three "families" of 
perturbations: 

1) The thermal operator e% and all the members of its sub-algebra (higher thermal operators). 



As shown in ||18|| , the Z^s PF theory perturbed by e\ flows into a massive regime, irrespec- 
tively of the sign of the coupling. Assuming that, as for the Z 2 case, due to the sign of the 
coupling Ai — A 2 in (|18D the theory is driven into a low temperature ordered phase, we have 
that vacuum expectation values (v.e.v.'s) of disorder operators fij, vanish for j ^ 2S mod{2S) 
as well as v.e.v.'s of the parafermionic fields (ipk^l) = 0, for 2k ^ 2S mod(2S). 
2) The family of the disorder operator (with the corresponding vertex operators) given by: 

[s] 

J2^ k :e^S Mz) - m) :+h.c. 

k=l 

where [S] means integer part of S. Note that p,2S corresponds to the identity operator and 
then for S integer, this family contains the single vertex operator 

: e ^(vW-^W) : +h .c. (19) 
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3) The family of the parafermionic fields: 



2S 



i k 



Y^Mzm(z) ■ e-JSaM*™*" : +h.c. 
fc=i 

Note that the system is invariant under the extended Z2S x Z2S transformation 

^ _> wPim-n)^ . w 2 P ( m -n )?/; ^t . 

\/27rm _ _ \f2ixn 

with w = exp(m/S) and m, n e Z. 

Since the parafermionic sector is massive, the effective theory for large scales can be 
obtained by integrating out these degrees of freedom. Since this massive sector is driven 
into the phase where (|20|) is unbroken, we can obtain the most general effective action for the 
remaining U(l) field, where only the vertex operators invariant under fl2"0|) will be allowed 
to appear. We get 



Z eff = Jv<j) exp (- J K s {d^f+ 



2S \ 
a s I cos(-=(p -£)) + ••• , (21) 



for S half-integer, and 



Z eff = Jv<f> exp (- j K s (d^) 2 + 
r S 

Ps J COB(-=(<f - 0))+ 



2S \ 
a s I cos(-={<p -$)) + ••• ), (22) 



2S 

for S integer. Here the dots simply mean higher powers of the perturbing vertex operators 
and Ks is an effective constant arising from the OPE of vertex and parafermionic operators 
in the process of integration of the massive degrees of freedom. We immediately notice that 
for integer S there is an extra vertex operator coming from (|19D which is not present for 
half-integer S, and as we will see, this difference between integer and half-integer effective 
actions is crucial. 

Using ([!]) and (H) we can write the continuum expression of the original spin operator 
S(x) as 

S(x) = J R + J L + {-l) x tr(a/2(<S>W + ^^)). (23) 

One way to see whether the system is gapped or not is to study the behavior of the spin- 
spin correlation function at large scales. Since our original SU(2) WZW model is perturbed, 
correlation functions of the fundamental field will contain supplementary operators coming 
from the OPE between c]^ 1 / 2 ) and the perturbing fields. With the help of the fusion rules 
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(14]) it is easy to see that, for example, the effective alternating z-component of the spin 
operator containing the scalar field will be given by: 

£ a kf , k : e^3 Mz) ~ m) : +h.c, (24) 

k<2S, k odd 

where only odd k fields appear in the sum. Let us consider now separately the case of 
half-integer and integer S: For S half-integer, the operator $( s < ) is present in (^41), and we 
can easily check that, (since fi2S corresponds to the identity), this operator is simply given 
by 

e V2S W V > + h.C. 



The other operators in the series contain parafermionic disorder operators whose correlators 
will decay exponentially to zero at large scales. Thus, considering only the Gaussian part of 
(pip, we can show that the spin correlation functions at large scales behave like: 

< S z (x)S z (y) >~ (-l)(*-v)\x-y\- 2SKa 
< S + (x)S_(y) >~ {-l)^\x - y\~ l K 2SK s) (25) 

The fact that the SU(2) symmetry is unbroken at all scales fixes then the value of Ks to be 

K s = 1/(25) (26) 



For this value of Kg one can show that the perturbing operator in ([21]) is marginally irrele- 
vant. 

We conclude then that the large scale behavior of half-integer spin chains is given by 
the level 1 SU(2) WZW model with logarithmic corrections as for the spin 1/2 chain. An 
interesting extension of this analysis is to add an arbitrary dimerization to the chain. In the 
context of the WZW approach this corresponds to the addition of the field tr (I^ 1 / 2 ) + h.c. oc 
fix : ezvW^ 2 -* : -\-h.c. in the perturbing terms ||. Reproducing the same arguments as 
above, we see that this term has the effect of making the operator cos(^=(y? — tp)) appear 
in the effective action ([H|). Using (|26|) , we see that now this operator is highly relevant and 
will so produce a gap in the excitation spectrum unless a fine tuning of the parameters is 
performed. 

Let us consider now integer spins S. Since the series (E3j) for the effective spin operator 
contains only half-integer spins j (odd fc's), all the operators in the series will contain non- 
trivial parafermionic operators. Then all the terms in the spin-spin correlation function 
will decay exponentially to zero with the distance indicating the presence of a gap in the 
excitation spectrum, thus confirming Haldane's conjecture. This is our main result. 

A possible modification consists in the addition of dimerization to the system (which 
again corresponds to the inclusion of the field tr as a perturbation). Then, also 

integer spin j terms (even fc's) will be generated (as indicated by the fusion rules (|T4])) 
in the effective spin operator and in particular the vertex operator flT9|) which is the only 
candidate for a power law decay of the spin-spin correlators. However, as for the (dimerized) 
half-integer spin case, the presence in ( f22|) of the (relevant) first vertex operator prevent the 
system to be massless (this situation is similar to the one encountered in [I9|, in the context 



of abelian bosonization), indicating again that a fine tuning of the dimerization parameter 
has to be performed to get a massless regime. 
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We have presented the continuum limit of the spin-S 1 antiferromagnetic Heisenberg chain 
in terms of a parafermion Conformal Field Theory interacting with a compact U(l) boson. 
After integrating out the massive parafermions we have shown that HAF chains with half- 
odd-integer spins have a massless spectrum while those with integer spin have a gap to all 
excitations, in complete accordance with Haldane's conjecture. 

It has been shown by using a mapping to the sigma model that a dimerized spin S chain 
should have 2S + 1 massless points in the dimerization parameter space ||. This result was 
shown to be valid for large S and a detailed treatment of dimerization within the present 
approach could help to extend this result for small values of the spin. 

We are grateful to F. Alcaraz, P. Dorey, V.S. Dotsenko, A. Lugo, V.A. Fateev, A. Ho- 
necker, M. Picco, V. Rittenberg for useful discussions. D.C.C. thanks CONICET and Fun- 
dacion Antorchas for financial support. C.v.R. is supported by CONICET, Argentina. 
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